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We investigate the dynamics of the Higgs field at the end of inflation in the minimal scenario 
consisting of an inflaton field coupled to the Standard Model only through the non-minimal gravita¬ 
tional coupling ^ of the Higgs field. Such a coupling is required by renormalisation of the Standard 
Model in curved space, and in the current scenario also by vacuum stability during high-scale in¬ 
flation. We find that for ^ > 1, rapidly changing spacetime curvature at the end of inflation leads 
to significant production of Higgs particles, potentially triggering a transition to a negative-energy 
Planck scale vacuum state and causing an immediate collapse of the Universe. 
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The Standard Model (SM) of particle physics can be 
consistently extrapolated to the Planck scale without 
any new physics, but the current measurements of the 
Higgs boson and top quark masses suggest that the cur¬ 
rent vacuum state of the Universe would then not be 
stable. This instability depends sensitively on the top 
mass TTit, which is subject to significant experimental and 
theoretical uncertainty 0 , but for the best fit values, 
the Higgs potential turns negative above the instability 
scale AI ~ lO^^GeV 2l|. This implies that the current 
vacuum would eventually decay into a negative-energy 
Planck scale true vacuum, but its lifetime exceeds the 
age of the universe by a wide margin [l|. 


Whether such a metastable universe could have sur¬ 
vived the cosmological evolution, especially inflation, has 
recently attracted significant interest [sj-y. In most of 
the simplest models of inflation, the Hubble rate dur¬ 
ing inflation is cormarable to the current upper bound 
H <9x 10^^ GeV [8|. It may therefore well be above the 
instability scale, in which case production of Higgs fluc¬ 
tuations could push the field over the potential barrier 
into the true Planck-scale vacuum Q. This instability 
problem is exacerbated by spacetime curvature induced 
running of the couplings, which makes the Higgs self- 
coupling negative even at low field values [3-0 • 


Notably, vacuum stability can still be maintained even 
during inflation without any new physics coupled to the 
SM fields [3], thanks to the Higgs-curvature coupling 
H. This coupling is inevitably generated by radia¬ 
tive corrections and when assuming the SM to be valid up 
to the Planck scale it is the only relevant new term when 
probing sub-Planckian scales. The current experimental 
constraints are extremely weak, |^| < 2.6 x 10^® Q. With 
a positive coupling, this term increases the height of the 
potential barrier between the vacua, thereby increasing 
the lifetime of the metastable vacuum. Vacuum stability 
is maintained for all inflationary scales compatible with 


the tensor bound (8| , provided the electroweak scale value 
of the running coupling lies above ^ew ^ 0.1 [3]. 

In this letter, we investigate the instability problem 
at the end of inflation, again assuming no new physics 
or higher-dimensional operators but taking the gravita¬ 
tional coupling ^ into account. In contrast with the Higgs 
inflation scenario 0, we assume that the Higgs is sub¬ 
dominant, and for simplicity we ignore a possible direct 
coupling to the inflation. We focus on the parametric 
region ^ = 1/6 in which case the Higgs is effec¬ 

tively massive during inflation and does not get displaced 
from the SM vacuum. When inflation ends, the curva¬ 
ture scalar R drops rapidly, reducing the height of the 
potential barrier. Eventually, a new barrier is generated 
by the thermalized degrees of freedom re-stabilising the 
vacuum. During the intermediate period rapid changes 
in the scalar curvature R can lead to non-adiabatic evo¬ 
lution, giving rise to significant excitations of the Higgs 
field, and potentially triggering the instability. Two ma¬ 
jor sources of excitations are a sudden drop in R when 
inflation ends iul, and parametric resonance from the os¬ 
cillating curvature term sourced by the inflaton 0[il. 

We can write the Standard Model Higgs Lagrangian as 


C = \D^H\^ - + ^R) \H\^ - X\H\^ 


( 1 ) 


where is the standard SU{2) x t/(l) covariant deriva¬ 
tive generalized to curved space, M the Higgs field mass 
parameter and ^ the nonminimal coupling. We will 
parametrize the complex Higgs field doublet H as 



hi + ih2 
hs + ihi 


( 2 ) 


Assuming the Friedmann-Robertson-Walker metric 
= dt^ — a?d-x? = a?'{drf' — dx^), and using the con¬ 
formal time coordinate rj, we write each of the four com- 





2 


ponents hi as 


Bogoliubov coefficients 


^^2TT)^a2 + 4 , (3) 

where k is the comoving momentum and “ 

(5(3)(k — = [®L“k'] = Ignoring interac¬ 

tions and the Minkowski space mass term , we obtain 
the mode equation 


Ok 


y/^out 


/Wir 



'\/^out 


/Wir 



For w ^ 1 /3, we can approximate Wi„/out = ™in/out in the 
long-wavelength limit k^ <C In that case the Bo¬ 

goliubov coefficients are independent of momentum and 
we find the occupation number 


/"(f) 


k“+(f 



fih) = 0, 


(4) 


rik = (in, Olfik^'^lOdn) = |/?kP 


(2 - Vl-3w)2 

8V1 - 3w 


where the prime denotes d/dri and R = 6a"/a^. 


In the matter-dominated case {w = 0) this reduces to 

Uk = 1/8. 


We can see that the mode has an effective curvature- 
induced mass term 

niLrv = 

If the frequency of the mode = k^ -|- satisfies 

the adiabaticity condition |w'^/a;^| ^ |a;"/w^| <C 1, the 
mode equation may be solved adiabatically via a WKB- 
type ansatz 


In the radiation-dominated case {w = 1/3), Wout = |k|. 
At long wavelengths, k^ ^ li'^’I ^ momentum- 

dependent occupation number 

IQ |2 

rik = |/3k| = —2|k|—’ 

where the star * refers to values at the time of the tran¬ 
sition and we have assumed ^ 1. 


fiv) = 


= dr]'uj 

\f2bJ 


( 6 ) 


The occupation number of a mode is an adiabatic invari¬ 
ant 1^, so if a mode starts in the vacuum and evolves 
adiabatically, it will stay unexcited. Only for a fast 
change in R do we expect to find sizeable excitations. 


For simplicity, we will first assume that the equation of 
state w = p/p changes discontinuously from w = Win = 
— 1 to a constant value w = Wont > 0 at the end of 
inflation. For constant w, the scale factor behaves as 
a = j and the scalar curvature expressed 

in terms of the Hubble rate H = d/a = a!/a} is i? = 3(1 — 
3w)H^. We assume that C ^ I; so that the adiabaticity 
condition is satisfied on both sides of the transition, but 
at the transition itself the curvature mass term drops 
instantaneously by a factor 


1 3Wout 


mt 


(7) 


This non-adiabatic change produces excitations of the 
Higgs field. 


Following [1^ we can write two sets of solutions: finiv) 
that reduces to the incoming adiabatic mode in the past 
and foutiv) that reduces to the outgoing mode in the 
future with = k^ and = k^ respec¬ 

tively. The two vacua can be connected with a Bogoli¬ 
ubov transformation finiv) = akfoutiv) + PkfLtiv) ,with 


In reality, the change in the equation of state is not 
instantaneous, and the slower the transition, the fewer 
modes are excited. To quantify this, we make use 
of Ref. [m (see also Ref. [1^), and model the time- 
dependence of the curvature mass term by 


riv) = 


+ ml,, 


mr„ - m- 


■tanhi^r;. (11) 


Here iz controls the speed of the transition, with v ^ oo 
and ^ 0 corresponding to infinitely fast and infinitely 

slow transitions, respectively. This case can also be 
solved analytically, and gives the occupation number 

_ sinh^ [tt (wout - Wjn) /(2tz)] 
sinh(7rwin/!z) sinh(7rWout/i^) 

From the various limits of (1121) follow a number of phys¬ 
ical implications: In the limit of large momentum, large 
mass or slow transition, i.e. k —> oo, ^ > oo or iz —> 0 the 
occupation number approaches zero exponentially. For a 
mode to receive excitations one must satisfy (assuming 

Win > Wout) 

|k| < min jmin, (zz/27r)2 - | . (13) 


The excitation of the Higgs field due to the non- 
adiabatic behaviour at the end of inflation leads to a 
rapid growth of the fluctuations of the field. To describe 
this we calculate the variance of the field in the ” in” vac- 
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uum state, at time rj after the transition 

(iihOl^LtlO.in) 

1 + 2|/3kP + 2ak/3kCOS ^ / drj'wj . 

(14) 

These fluctuations are uncorrelated on superhorizon 
scales, and therefore they correspond to a stochastic 
background as opposed to the coherent Higgs condensate 
generated during inflation 0- This will lead to vacuum 
decay, if the field value averaged over a sufficiently large 
volume, which we conservatively consider to be the whole 
Hubble volume, exceeds the position of the potential bar¬ 
rier. Therefore we only include modes up to A = aH. We 
also neglect the vacuum part and coherently oscillating 
terms in da. 


d3|k| 


2w(27r)3 


{h^)A 


d|k|k^|/3kP 


(15) 


Substituting Eqs. ([9]) and (flOll . we get for ^ » 1 



itive and negative values. As a result, the Higgs field 
grows exponentially in a way that can be interpreted ei¬ 
ther as tachyonic or resonant growth. 

Reheatin g vi a the curvature coupling ^ was first stud¬ 
ied in Ref. (1^ where it was named ’geometric reheat¬ 
ing’ (see also [13 )■ A similar tachyonic resonance effect 
present in reheating models with trilinear interactions 
was studied in 17| . Assuming that the inflaton potential 
can be well approximated by the quadratic term 




(18) 


we can write the solution as = ^cos{mt), where 
$ ft! V6HMpi/ m is a slowly varying amplitude and 
H = 2/(3t) as given by the equation of state for mat¬ 
ter. Substituting this background solution into Eq. (El) 
we obtain the Mathieu form of the mode equation ([4]): 


dz^ 



2qcos{2z) 


f{z) = 0, 


z = mt, (19) 


Ak = 


k2 






2A/pV 


9 = 


3$^ 

4A^ 




Following the analysis of [17| we get for the occupation 
number after the first oscillation 


From (El) and © we see that for w ^ 1/3 the super¬ 
horizon variance is suppressed by This is however 

not the case for conformal equation of state w = 1/3 due 
to inverse scaling whereby a large variance may be 
generated for ^ ^ 1 for a sufficiently rapid drop in mcurv 

As an example, we assume a model where inflation 
ends within a time scale ^ (400iJ)“^ with an effec¬ 
tive reheating equation of state w = 1/3. Using (fTOl) with 
^ > 500 we get for the magnitude of the generated vari¬ 
ance: Therefore we may conclude that if 

inflation ends abruptly to a state with rcout ^ 1/3, a vari¬ 
ance larger than A/ can be generated for sufficiently high 
inflationary scale H > Aj, such that the fatal transition 
to the negative energy vacuum is likely triggered. 

Above we assumed that the equation of state changes 
monotonically to its final value when inflation ends. How¬ 
ever, this is not what happens in most models of inflation. 
Instead, the inflaton field oscillates coherently about its 
minimum, and this leads to a much stronger effect. For 
an inflationary model with a single inflaton field tp in po¬ 
tential V{(j)), the scalar curvature R is given by 


= Xk= / Hkdzftv^^, (20) 

Jaz 

where with oj'^ being the term in the square 

brackets in m and Az covers the time when < 0 
during first oscillation. In the last form of (1^01) we have 
assumed that |k| < aH and that the amplitude $ is 
roughly constant (a ^ 1) during the oscillation. It is 
clear from (12(11) that for large ^ the occupation number 
may become very large already during the first oscilla¬ 
tion. Using (I^Hl) we get an estimate for the superhorizon 
variance after the first oscillation (|/3kP = 

^ d|k|k2|/3kP ^ /gy 2exp{v/?^} 
27r2wout V27ry 

( 21 ) 

which is accurate for ^ > 1/6-I-3/8 ~ 0.5. Below 
this region only superhorizon modes with k^ < = 

— 3/8) receive tachyonic amplihcation. In this 
case, an accurate approximation for the superhorizon 
variance can be obtained by an expansion near the 
threshold point ^ = 3/8 




It is clear from this expression that i?, and therefore also 
the curvature mass term oscillates between pos- 


( 22 ) 

The exponential growth of the particle number and 
the variance is constrained by backreaction which even- 
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tually shuts off the tachyonic resonance and makes the 
dynamics non-linear. We consider two backreaction ef¬ 
fects: Higgs self-interaction and gravity. 

When considering the Higgs self-interaction, we adopt 
the convention that the RG-running parameters are eval¬ 
uated at the scale ^ = H, which corresponds to the op¬ 
timal choice in terms of RG-improved effective potential 
0 in case H is the highest scale of the problem [2^. The 
RG-running self-coupling A(/i) becomes negative at the 
scale of instability A/ and therefore we may approximate 
= H): 

X{H) ~ Ao sign(A 7 — H ), with Aq ~ 0.01. (23) 


Hence, the backreaction from self-interactions is limited 
to the region H < Aj, while for H > Aj the negative 
self-coupling enhances the tachyonic resonance and ac¬ 
celerates particle production. However, here we will not 
consider this enhancement but use the conservative esti¬ 
mate of m- 

For H < Ai, the positive Higgs self-interaction gener¬ 
ates an effective mass 


"left = 'mlurv 


6Ao(h2). 


(24) 


The resonance is shut off once (see also [1^ [l^) Ak + 
6A > 2\q\ with 5A = 6X{h^)/m^. From this and 
Eqs. (I23II we then obtain a condition for maximal vari¬ 
ance in the region H < A/: (h^) < ^iJ^/Ao, including 
contributions from all tachyonic modes |k| < fccut- The 
corresponding superhorizon variance with |k| < aH is re¬ 
lated to the total variance by {h'^)aH ^ 
and is hence constrained by 






aH 


6AoV3? 


for 


H <A 


1 1 


(25) 


where we have used R ~ for ic r; 0 assumed to be 
valid after oscillations (or in average during oscillations). 

For gravity backreaction we need to compare the en¬ 
ergy density contained in produced Higgs particles to to¬ 
tal energy density p ~ 3i/^Mp[. In Gaussian (Hartree) 
approximation with (h^) = 3{h^)^ we get for the Higgs 
energy density an estimate 


PHiggs 2^H^{h^) + 6X{hY , (26) 


whereby we find the following constraints for the super¬ 
horizon variance for the Higgs energy density to remain 
sub-dominant during oscillations |24j : 


{h^)aH ^ 


2Mp2 

ll(3C)5/2 


and (h^}aH < 


HAL 


pi 


2(2Ao)i/2(303/2 ■ 
(27) 

We estimate that the probability for a potentially 
catastrophic transition to the negative energy vacuum is 



FIG. 1: The estimated instability region (shaded, blue) 
where Ah > Aj after the first oscillation, given by 
Eas. (l21l22l25l27l28ll . The dashed lines correspond to Ah > 
lOAj and Ah > 10^Aj. For the parameters we have used 
4? = O.SMpi and Aj = 10“^Mpi. 


significant once the superhorizon variance has exceeded 
the instability scale A/ ^ 10^^ GeV of the Higgs effective 
potential: 


Ah = \J{\H\^)aH = \J‘^{h‘^)aH > Aj . (28) 

We show the corresponding region after the first oscilla¬ 
tion as blue in Figure [T] On the left, at ^ < 10^, the 
increase of the variance is limited by the weakness of 
the resonance (l2T|l . although the non-shaded area above 
the dashed blue line would be mostly filled by the subse¬ 
quent 3-4 oscillations. Further to the left, the effect be¬ 
comes even weaker and disappears completely at ^ = 3/8. 
The bottom half of the plot is excluded by backreaction 
from the self-interactions. On the right, the backreac¬ 
tion from the gravity effects makes the linear approxima¬ 
tion invalid, and more sophisticated methods are needed. 
Backreaction from the other SM degrees of freedom is not 
likely to change the qualitative picture as resonant ampli¬ 
fication takes at least a few inflaton oscillations. A direct 
Higgs-inflaton coupling can enhance or weaken the effect 
depending on the model but would generally not elimi¬ 
nate it as long as the Higgs has a tachyonic phase during 
the hrst inflaton oscillation. 

In conclusion, we have shown that for a sufficiently 
high inflationary scale H > Aj, changing spacetime cur¬ 
vature can trigger a fatal transition to a negative en¬ 
ergy vacuum, for a wide range of the curvature coupling 
1 <^ < 10^. This conclusion applies to both rapid re¬ 
heating and oscillating inflaton. Gombining this with the 
lower bound ^ew ^ O.I imposed by stability during in¬ 
flation , we find that the minimal scenario of Standard 
Model in a background of conventional high-scale infla¬ 
tion constrains the coupling ^ to be close to its conformal 
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value ^ = 1/6. 
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